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Abstract 

Solutions of RG equations for f3(a) and a{Q) are found in the class 
of meromorphic functions satisfying asymptotic conditions at large Q 
(resp. small a), and analyticity properties in the Q 2 plane. The 
resulting (Xr(Q) is finite in the Euclidean Q 2 region and agrees well 
at Q > 1 GeV with the WS a s {Q). 



1 Introduction 

QCD is known to simplify at large N c : i) the perturbation theory dictates 
that only planar diagrams survive in the leading order pQ; ii) assuming that 
confinement exists also in the limit N c — > oo, the spectrum of mesons and 
glueballs consists of bound-state poles, and the decay width vanishes at large 

In the 1 + 1 QCD this property was proved both analytically and numer- 
ically [3] (for a review see 0]) since confinement in this case is induced by 
the perturbative gluon exchange. 

On the lattice numerical calculations have confirmed that the large N t 
limit is achieved with few percent accuracy already at N c = 3 for several 
tested observables [3]. It is thus likely that the large N c QCD is a good first 
approximation for the realistic N c = 3 case, which has an accuracy of 10% 
or better. 
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From theoretical point of view the large N c limit of QCD is very useful 
since it drastically simplifies the analytic structure of amplitudes, e.g. one 
has sums of simple poles in the two-point function, and for the four-point 
function one expects formulas of the Veneziano type. 

Thus one expects the large N c amplitudes as functions of external mo- 
menta to be meromorphic. On the other hand the perturbation series yields 
the typical logarithmic dependencies and the RG equation prescribes for 
a s (Q) the structure with unphysical poles and cuts which are incompatible 
with unitarity and analyticity. E.g. the one- loop expression for a s (Q) has 
the form 

47T 11 2 

" s(g) = /3olngVA 2 ' ^ = ~3 N ° ~ 3^ (1) 
with the Landau ghost pole at Q 2 = A 2 and the cut with the branch point 
Q 2 = 0. This property being true for all N c , strikingly violates the expected 
meromorphic structure of amplitudes for iV c — > oo. 

One might argue that a s (Q) by itself is not yet the physical amplitude, 
and in the latter the unphysical features of a s may be somehow compensated. 
This is however contradicted by examples of amplitudes, e.g. for the process 
e + e _ — > hadrons, where enters directly into the hadronic ratio R(s). 

R(s)=R parton + ^ + o(^) 2 . (2) 

Moreover one can define the "effective coupling" a R (Q) for the process R, 
which enters the physical amplitude Q R directly 

n R = {l { ° ] +u R a R (Q) (3) 
and this "process dependent" a R satisfies the standard RG equation 

da ^ = 8( R )(a) = pi ° R)a2 - ^ )a3 + 0(a 4 ) (4) 

where /3q R ^ and (3^ are standard scheme-independent coefficients, while 
P^\n > 1 depend on the process. One can easily see, that the solution 
of (@J) has the dominant asymptotic term of the form (0) with the unphysical 
features discussed above. 

One might still argue that perturbative series and a s (Q) itself should be 
considered only in the asymptotic regime when Q is large and therefore the 
logarithmic singularities and Landau ghost pole are far away. 



2 



However the analytic structure of Riemann surfaces with cuts is different 
everywhere from that of meromorphic function and this argument of asymp- 
totics can be turned around to imply that the logarithmic asymptotics of 
a s (Q) is an asymptotic approximation of the true analytic function which is 
meromorphic at large N c and is in agreement with all expected properties of 
large N c physical amplitudes. 

It is the purpose of the present paper to exemplify the solutions of the 
RG equations which have the desired meromorphic properties. These solu- 
tions will have the standard logarithmic asymptotics in good numerical and 
analytic similarity with the standard perturbation theory. 

Moreover we show that this meromorphic-logarithmic duality has deeper 
roots in the quark-hadron duality and provide explicit example of this con- 
nection. 

The paper is organized as follows. In section 2 the RG equation for a s (Q) 
is considered and the general form of solution is written down, containing an 
arbitrary function with known asymptotics. In section 3 the meromorphic 
solution is suggested and its properties are studies. In section 4 connection 
with the standard perturbation theory is investigated. The concluding sec- 
tion is devoted to general discussion of the meromorphic - logarithmic duality 
- in QCD and QED. 



2 General solutions of RG equations 

In this section we follow the line of reasoning which was given in [7j. We 
write the RG equation as in (HJ) suppressing the subscript R everywhere and 
express (3(a) through another unknown function <p 

da(n) fa a 2 

P(a), p(a) = - — j- r . (5) 



Here (p'(x) means the derivative of function ip(x) in the argument x, and 
fi is the RG scale, which will be later traded as usual for external parameters 
(momenta) of the given process Pf since a can depend only on the ratio 

/?//'• 

In terms of the function (p the solution of the RG equation for a(fi) can 
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be immediately written 

4:71 



a: 



(6) 



/3 Q [ln^ + X + ^(l);" 

Here \ is an arbitrary function of Pi such that In fi 2 + \ is some function of 
P 2 I fj, 2 , P 2 / P 2 the form of which depends on the process and will be found 
below for concrete examples. 

The Standard Perturbation Theory (SPT) which is assumed to be valid 
at large Q 2 , provides some limitations on the properties of functions \ and 
ip, which one must impose. Namely, the SPT Taylor expansion of /3(a) is 
known for the first four terms 

a( \ ft 2 ft 3 ft 4 ft 5 . 6n fn\ 

a (a = a -a -a — - — —a+(J{a). 7 

HK ' 2tt 4tt 2 64tt 3 (4tt) 4 v ; v ' 

Here /3 , ft are scheme- independent and for rif = are equal to (3q = \-N c — 
|n/ and /3i = ^-N 2 — \N c rif — N c nf, while /3 2 and /3 3 have been calculated in 
the MS scheme. Comparing and ((HJ) one obtains the following condition 
on the asymptotic behaviour of <f(x) at large x(x = -): 

^)Uoo = lnx + oQ. (8) 

The first term on the r.h.s. of (JBJ) reproduces the scheme-independent co- 
efficients ft, ft in the expansion while the term O(-) contributes to 
the higher order coefficients ft, ft, ... Conditions on the function x{P 2 ) are 
more subtle and in general depend on the process in question. Here one 
should distinguish in Pf the external parameters which can be made arbi- 
trarily large, as the momentum Q 2 in the two-point function H(Q 2 ) or in the 
formfactor, and other renorm-invariant parameters, which define the scale of 
confinement, e.g. the string tension a, or the RG invariant gluonic conden- 
sate ^^(F Miy (0)F /iI/ (0)). In the framework of the Background Perturbation 
Theory (BPT) jHJ E] both appear as vacuum expectation values of operators 
made of RG invariant combinations gF^, where refers to the background 
field. Here we do not use BPT, but only consider all RG invariant parame- 
ters like a on the same ground as the true external parameters - the external 
momenta like Q 2 . In what follows we shall use the generic mass parameter 
m 2 = Att<7 as the confinement scale and disregard for simplicity other vac- 
uum field characteristics, like the gluonic correlation length A, which can be 
computed in principle and in practice PjH 3 through m 2 . 
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As a result keeping only one external momentum Q 2 , x can be written 
as x = x(Q 2 / m2 )- Then the large Q 2 behaviour of a(Q 2 ), 



x 47r ^/lnln0 2 \ 

dictates the following behaviour of x+ m A* 2 ? where we go over to the standard 
Aqcd parametrization instead of \i parametrization 

m 2 O 2 
( X + ln/i 2 )=ln — + *(Q 2 ,m 2 ), s In ^. (10) 

3 Meromorphic realization of RG solutions 

As was discussed in Introduction, at large N c the analytic structure of phys- 
ical amplitudes and of effective charge a(Q) is simplified and reduced to the 
meromorphic functions with simple poles at Q 2 = —M 2 , corresponding to 
the bound states of quarks and gluons - mesons, glueballs and hybrids. 

Therefore the function ^(Q 2 ,m 2 ) in (}2T)j) can be represented in the spec- 
tral form 

nQ\rn 2 ) = -±-^— + const (11) 

n=0 V I 1V1 n 

where the coefficients C-yi £1X6 independent of Q 2 . The spectrum M 2 was 
found for large N c for mesons, glueballs and hybrids in the limit of no mixing 
between them jTJ] . In the large N c limit the meson-glueball mixing vanishes, 
while the meson- hybrid mixing is 0(a) jT5j, and following JT] we write the 
lowest order in N c and a-independent spectrum, for glueballs and mesons as 

M^ 2 = m 2 (n + L/2) + M 2 + 0(1/ n) (12) 

where m 2 = Ana, and for glueballs a(glue) = ^a(quark). In what follows 
we shall be interested only in the asymptotic part of the spectrum at large n 
and L = 0,1, (the corrections to (JT2J) for n = 1, 2 come from spin splittings 
and is of the order of 20%). In general one can rewrite the effective coupling 
for the given process R with external parameters denoted as Q 2 as 

An 

a R = — f ; : : — — ■ ( 13 ) 



0o 



ln^ + *«(Q 2 ,m|) + f^(i) 



5 



Here m R , ty R and ip R depend on the process and ^ R contains the poles of 
the spectrum in the lowest (one loop) approximation. Conditions on ^/ R , ip R , 
Eqs. (fTUj) and (JBJ) respectively, impose restrictions on coefficients c n , in 
the spectral representation (jllj) . As an example one can choose ^ R (Q 2 , m 2 R ) 
in the following form 

, (Q 2 + (M^) 2 \ s T'(x) 

\ m K / r(ar) 

which implies that the spectrum responsible for the one-loop RG evolution 
of a R (Q) is 

M 2 RN (1 - loop) = m\n + (M R 0) ) 2 . (15) 
The asymptotics of ip(x) at large x is 

^(x)=\nx- — -J2^k ( 16 ) 

where B 2 k are Bernoulli numbers behaving at large k as B 2 k ~ 2(— ) 1 (2A;)!/(27r) . 

Therefore the asymptotic condition ([TO]) is satisfied by the choice (JHJ). 
Note that one can also satisfy this condition by the realistic spectrum which 
asymptotically has the form (|TH|) but differs from it for the first Nq terms, 
provided the coefficients c n tend to constant for large n. 

At this point it is important to define the physical system which pro- 
vides the one-loop spectrum ([15)1 . To this end one can compare the one-loop 
expansion in SPT, a s (l - loop) = a< > - § \n(Q 2 / /j, 2 )(a^) 2 + ... with the 
equivalent expansion of a R (JT3J): 

a H ("l - loop") = c4 0) - £^(g 2 ,m|)(ag ) ) 2 + ... 

where a R = — One can see that ^ R (Q 2 , m R ) plays the role of gluon 

loop, or better, two-gluon intermediate states imbedded in the framework of 
the process R. In particular, when R is the e + e~ annihilation into hadrons, 
then *5/ R is responsible for bound states of quark, antiquark and two gluons, 
i.e. in the large N c limit this is the two-gluon hybrid state. 
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4 Meromorphic solution and Standard Per- 
turbation Theory 

In previous section it was argued that the general solution of RG equa- 
tion @ which correctly reproduces the one-loop result of SPT and scheme- 
independent coefficients /3 , /?i in the expansion of /3(a) is given by 

47T 

a R= f o 7^ ,..rn^N : — ■ ( 17 ) 
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where <Pr(x) — > Inx, x — > oo. The analytic properties of an(Q 2 ) depend 
on the form of tpn (J^\ . Neglecting the latter (one-loop approximation) one 
obtains for an the meromorphic function of Q 2 with simple poles at the 
negative values of Q 2 (for m 2 R > A 2 ). If one adopts (pn{x) = Inx, then one can 
reproduce the two-loop result of SPT, but then ccr acquires the logarithmic 
branch points, in the Q 2 plane which are unacceptable for N c — > oo. To 
keep the correct meromorphic properties of atR one can choose lpr as a 
meromorphic function ot its argument. A particular choice was made in [7], 
namely 



where ip(x) is the same as in (|14j) with the asymptotics which satisfies 
asymptotic condition on (Pr(x) given above. The resulting expression for 
0(a) is 

fl^-l x-a^ + AV (19) 

One can persuade oneself that for A > A = 1.255, 0(a) is analytic for 
a > 0, and has infinite number of zeros for negative a condensing on the 
negative side of the point a = 0. For positive a the function (3(a) is always 
negative and monotonically decreasing, and thus has no IR fixed point. Also 
the coefficients 02,03 computed from (|19|) are several times smaller than 
computed in MS scheme. 

To compare or in Eqs. (|17jl. (j!8)l with the loop expansion of SPT one 
can consider Q 2 large enough so that one can use the asymptotics (fTKJ) for 
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ip{x) and keep the leading (logarithmic) term. With the notation 



L = ln 



q + m 



(0)\2 



a R — 



Air 



A 2 ' - fcV 
one obtains the expansion in 1/L and which looks like 



(0) i 2 fa In L API 
a R = a> R ' \\ ^ — — + 



PI L ' Pil? 



In/.--) +b 



+ 



InL 



Here6=-f (A-i 



|, and for A 



(20) 



(21) 



1.5 one has b 



•1.436, while in the 



MS b = 0.26 for rif = 0. Thus the difference with SPT occurs only in the 
3-loop result and is not large. Typically at Q 2 = 3 GeV 2 the exact expression 
(jnj) yields a R = 0.26 while the 2-loop WS result is a R = 0.256(ra / = 0). 



5 Conclusions 

We conclude this paper with the discussion of the OPE for the case of e + e~ 
annihilation which will stress the new features of the a. R (Q 2 ) behaviour (|13|). 
Expanding the polarization function H(Q) in powers of a s one has 

n(Q) = n(°)(Q) +« s n«(Q) + a 2 Il^(Q) + ... (22) 

where U {0) {Q) can be computed as the spectral sum |T3] as follows 

n(o) (Q) _ j_ v ^ - f g!±ggn (23) 

and similarly for the hybrid sums 

This should be compared with the standard OPE expansion ^3] 

where m g = m "+ m rf . As was discussed in US], the leading asymptotic 
term in ()23)1 reproduces the partonic answer ^^Tn^-, supporting the idea 
of Quark-Hadron Duality (QHD). 
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The next asymptotic terms of (J23|) produce the OPE terms 0(^-^j, 
O \ Tjz)i etc. One usually imposes the QHD requirement, trying to reproduce 
the OPE coefficients in and in particular the absence of O (rjs) term 

for m q = ^31- Here we notice, that the form (J 17)) produces additional 
power terms of somewhat different structure. Namely, taking the subleading 
asymptotic term in (|T6*j) for an one has 

(°) , m2 ^° ( (OK 2 , / ot r\ 

aR = aR + 2(Q* + m^ {aR) + - (25) 

Assuming that UP) (Q) has the same asymptotics as Il( ) ( to reproduce the 
term ^ In % in (j2U)) one obtains a new power term 

AI1(Q 2 ) = - A— ^ (26) 
w ' 12n 2 Ti 2Q 2 v ; 

which is O { qz^qi ) anc ^ nas a negative sign. 

Thus the meromorphic a^, Eq. (|17|) produces the power-like terms of 
new kind. Numerically these terms are smaller than those originating from 
asymptotic expansion of ()23|). and phenomenologically the appearance of 
0(1/Q 2 ) terms is not forbidden by experiment |15j . 

As a summary, a new type of solutions of RG equations is suggested 
where an(Q) is a meromorphic function of Q with poles corresponding to 
the physical poles of two-gluon hybrid states. The resulting an(Q) is finite 
for Q 2 > and agrees well with ajjs(Q) for Q > 1 GeV and is phenomeno- 
logicaly acceptable for all positive Q 2 . 

The author is grateful to I.Ya.Arefieva, A.M.Badalian, A.L.Kataev, N. V.Kras- 
nikov, A.A.Pivovarov, D.V.Shirkov, for useful remarks and discussions. The 
work is supported by the Federal Program of the Russian Ministry of Indus- 
try, Science, and Technology No 40.052.1.1.1112. 
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